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1 Abstract 


Actin and myosin assemble into a thin layer of a highly dynamic network underneath the membrane of eukaryotic 
cells. This network generates the forces that drive cell- and tissue-scale morphogenetic processes. The effective 
material properties of this active network determine large-scale deformations and other morphogenetic events. For 
example, the characteristic time of stress relaxation (the Maxwell time tm) in the actomyosin sets the time scale of 
large-scale deformation of the cortex. Similarly, the characteristic length of stress propagation (the hydrodynamic 
length A) sets the length scale of slow deformations, and a large hydrodynamic length is a prerequisite for long- 
ranged cortical flows. Here we introduce a method to determine physical parameters of the actomyosin cortical 
layer in vivo. For this we investigate the relaxation dynamics of the cortex in response to laser ablation in the 
one-cell-stage C. elegans embryo and in the gastrulating zebrafish embryo. These responses can be interpreted 
using a coarse grained physical description of the cortex in terms of a two dimensional thin film of an active 
viscoelastic gel. To determine the Maxwell time tm, the hydrodynamic length A and the ratio of active stress 
and per-area friction 7 , we evaluated the response to laser ablation in two different ways: by quantifying flow and 
density fields as a function of space and time, and by determining the time evolution of the shape of the ablated 
region. Importantly, both methods provide best fit physical parameters that are in close agreement with each 
other and that are similar to previous estimates in the two systems. Our method provides an accurate and robust 
means for measuring physical parameters of the actomyosin cortical layer. It can be useful for investigations of 
actomyosin mechanics at the cellular-scale, but also for providing insights in the active mechanics processes that 
govern tissue-scale morphogenesis. 

2 Introduction 

Cells need to adopt their shape in order to drive tissue scale morphogenetic processes lUElEllll. A classical 
example is convergent extension, where cells reshape and intercalate in order to allow an epithelium to shrink in 
one direction, while expanding in the other 15]. The actomyosin cortex also endows cells with the ability to reshape 
themselves la. This thin layer beneath the membrane largely consists of cross-linked actin filaments and non¬ 
muscle myosin motor proteins. Importantly, this thin structure can generate active stresses and contract EElil. 
Active stresses emerge from the force-generation of myosin motors interacting with actin filaments, fuelled by 
ATP hydrolysis ifTOl . Such active molecular processes lead to the buildup of mechanical stress (cr®^) on larger 
scales im. To gain a physical insight into the stresses and the force balances that govern large-scale deformation 
and flows of the cell cortex, coarse grained continuum mechanical descriptions have played an important role 
iiiiiiii. Furthermore, cortical laser ablation (COLA) has emerged as a useful tool for investigating forces and 
stresses in the cortical layer |l3][8l|3. The aim of this manuscript is to extend the analysis of the response of the 
actomyosin cortex to COLA by use of thin film active viscoelastic gel theory iniiTiiiTS], in order to determine 
physical parameters that characterize the actomyosin cortical layer. 
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One biological example where a continuum mechanics description in the framework of active gels has been 
useful is epiboly in zebrafish gastrulation ||8]. Here, a ring of actomyosin cortex forms on the surface of the 
yolk cell ifTbl . This ring contracts to generate circumferential stresses, but also to drive flow of actomyosin into 
the ring. The mechanical stresses that are generated in this process play an important role to pull the connected 
enveloping layer of cells (EVL) of the blastomere from the animal pole of the embryo towards the vegetal one. The 
stress distribution and the emerging large-scale actomyosin flow velocity field can be understood using thin active 
gel theory. Cortical tension can also be investigated in experimental terms by cortical laser ablation (COLA). 
Anisotropies of recoil velocities can be related to tension anisotropies resulting from flow, and viscous tension and 
the emerging flow fields depend on the large-scale physical parameters of the actomyosin cortical layer. Therefore, 
in order to understand the mechanics that underlie flow and deformation of the actomysoin cortical layer it is key 
to determine its large scale physical properties, which has been difficult to achieve. 

A second example that highlights the role of actomyosin in morphogenesis is the polarizing one-cell stage 
C. elegans embryo ElEl. Here, gradients of active tension generate cortical flows that lead to the establishment 
of anterior-posterior cell polarity, which in turn is key for the subsequent asymmetric cell division imilQl. No¬ 
tably, cortical flow in this system can be well described by a thin film active gel theory 0. Also here, COLA has 
permitted the characterization of tension profiles 0. However, an estimation of other physical parameters that 
characterize the actomyosin cortex, such as the effective 2D viscosity p, the friction coefficient 7 with respect to 
the membrane and the cytosol have remained elusive. Furthermore, the cortex behaves as an elastic solid on short 
times, while it is essentially viscous on long times ED. Thus, a dynamic description needs to take into account 
that the actomyosin cortical layer is viscoelastic. A simple model for viscoelastic behaviors which is limited to 
the slowest relaxation processes is the Maxwell model which can be incorporated in the theory of active gels 
iniiiiiTa. Viscoelastic behavior of the cortex stems from the continuous remodeling of the cortical network 
by the turnover of actin filaments and other cortical constituents ll22l \23\ l2^ . This turnover occurs under the 
influence of regulatory and signaling molecules Col . Because of turnover, the strain energy stored in the cortex 
elastic stress relaxes and the corresponding strain energy is dissipated ll^l26l . The corresponding time scale of 
stress relaxation (the Maxwell time tm) is an important physical quantity that determines cortex behavior. Me¬ 
chanical perturbations that persist on times large compared to tm will lead to viscous deformations and flows 0. 
If, however a mechanical perturbation persist only for times that are shorter than tm, strain energy will be stored 
and the the cortical layer will respond elastically. To conclude, the time scale of cortical remodelling determines 
the characteristic timescale of stress relaxation tm, which in turn governs the cortical response to mechanical 
perturbations. 

An second important parameter characterizing the coarse-grained spatiotemporal dynamics of the actomyosin 
cortex is its hydrodynamic length, A = This length determines the range of stress propagation and sets the 

correlation length of the cortical flow field. Hence, a large hydrodynamic length A leads to long ranged cortical 
flows, which, for the case of C. elegans, is important for polarization GOl . In this manuscript we sought to 
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determine key coarse grained physical parameters of the cell cortex , by use of COLA experiments and by use of 
a theoretical description of the cortical layer in terms of a two dimensional active viscoelastic gel. We apply our 
method to two model systems: the actomyosin ring that drives zebrafish epiboly, and the actomyosin cortical layer 
in the single cell embryo of C. elegans that drives cell polarization. We measure distinct physical properties of the 
cortex in the two systems, thus demonstrating a broad applicability of our method. 

3 Materials and Methods 

3.1 C. elegans strain and sample preparation 

To image NMY-2 in single cell embryos, we used transgenic line LP133 (nmy-2(cp8[NMY-2::GFP + unc-119(+)])!; 
unc-119 (ed3) III). C. elegans maintenance and handling was as previously described lIZTl . We cultured the line 
at 20 C° and shifted temperature up to 24 C° 24 hours prior to microscope imaging and COLA. Embryos were 
dissected in M9 buffer (0.1 M NaCl and 4 % sucrose) and mounted onto the agar pads (2 % agarose in water) to 
squish the embryos gently. All experiments were performed at 23 - 24 C°. 

3.2 Zebrafish transgenic lines and sample preparation 

To visualize NMY-2 in the yolk syntial layer (YSL) of zebrafish embryos throughout epiboly, we used transgenic 
line Tg(actb2:myll2.1-EGFP) ll28l . Maintenance of the fish line as well as embryo collection were conducted 
as previously described ||29]| . Zebrafish embryos were incubated at 25-31°C in E3 medium and staged according 
to morphological criteria ||30l. For imaging and the performance of COLA embryos were mounted in 1 % low 
melting point agarose (Invitrogen) inside E3 medium on a glass bottom petri dish (MatTek). 

3.3 Imaging and COLA 

Zebrafish embryos and C. elegans embryos were imaged and ablated using modified versions of previously de¬ 
scribed spinning disk laser ablation systems 10 El. In brief, the spinning disc system (Andor Revolution Imaging 
System, Yokogawa CSU-Xl) was assembled onto a the Axio Observer Z1 inverted microscope (Zeiss) equipped 
with an 63x water immersion objective. Fluorescent images were acquired by an EMCCD (Andor iXon) at the 
specified time intervals (for C. elegans: 1 sec, for zebrafish 0.5 sec). The pulsed 355nm-UV laser (Powerchip, 
Teem Photonics) with a repetition rate of IkHz was coupled into the Axio Observer and steered for point-wise 
ablation by galvanometric mirrors (Lightning DS, Cambridge Technology). A custom-built Lab View program 
integrated all devices for simultaneous COLA and imaging. To see the cortical response in C. elegans embryos, 
we applied 20 pulses for each points spaced every 0.5 /rm along a 10 pm line. The cut lines were chosen to 
be parallel with the long axis of the embryos, the future AP axis. For zebrafish embryos we applied 25 pulses 
per point spaced at 0.5 pm along a 20 pm line. The cut line was placed within the YSL actomyosin ring at a 
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distance of 20 /j,to from the EVL margin and parallel to that margin. The intensity of the UV laser was adjusted 
to achieve successful COLA in zebrafish embryos. Successful COLA was characterized by the visible opening 
of the cortex in response to the cut, subsequent recovery of the actomyosin cortex within the cut opening and no 
‘wound healing’ response, as previously described iiisi. 


3.4 Comparison to theory 

We convert the NMY-2 fluorescence intensity I{x, y, t) to a scalar height held h{x, y, t) according to 


h{x,y,t) 


- Ip 

-^max(^) ^0 


( 1 ) 


Here, Iq denotes the ‘background’ signal obtained from the average intensity within the box located at the center 
of the cut opening in the hrst post-cut frame, /niax(f) is the maximum intensity in each recorded image, and x and 
y denoting spatial cartesian coordinates within the plane of the cortical layer and t denoting time. Note that this 
conversion assumes that myosin intensity is proportional to cortex height. 

To determine the best ht non-dimensionalized model parameters ai, 02 , 0:3 and the characteristic time Ta 
(see Appendix for the detail), we performed the least square htting for the COLA responses in experimental 
observation by theoretical one. We evaluated the spatial velocity prohle at a time just following the COLA, 
temporal evolution of the width of the cut opening boundary, and the recovery timecourse (see Lig. [^. Litting 
was performed iteratively by using the Nelder-Mead method ED. Lirst we obtained the best fit value of Ta with 
the arbitrary values of a's by htting the recovery timecourse. We then nondimensionalized the experimental time 
course with determined Ta value and obtained the best ht values of a's through htting the spatial velocity prohle 
and the time evolution of the cut opening boundary. In next iteration step, we determined the best ht value of 
with a's values obtained in the previous iteration step. We stopped the iteration after the convergence of the htting 
parameters. 

In parallel to this, we tested to hnd the best ht parameter values by comparing the temporal development of the 
cut boundary between the experimental data and theoretical prediction. We determined the cut opening boundary 
by automatically detecting the cut opening region. The edge points were detected by using active contour method, 
a built-in function in Matlab (Mathworks). Then we obtained the best ht parameter values iteratively by using the 
Nelder-Mead method ED- In each iterations, we hrst determined the best ht value of Ta by htting the recovery 
curve with the arbitrary values of a's, to set the time scale for the non-dimensionalization. Then we computed the 
cut response numerically with the given parameter values. We next compared the edge points of the cut opening 
region between the experimental observation, {xf{n),yf{n)), and theoretical prediction, (a:*(n), j/*(n)), at the 
given time point, n. We computed the pairwise distances between (a;®(n), j/®(n)) and (a;‘(n), ?/j(n)), dij{n), and 
then calculated the minimum distances, with respect to each j, df^{n) = mmj[dij{n)]. Then we took sum with 
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respect to i, and average for the frames in analysis, to get the distance measure, 


^ = EE <'(")■ ® 

n i 

where N represents the number of frames in analysis. We minimized D to obtain the best fit parameter values of 
a’s with the given value of Ta. The iteration is stopped after the convergence of the parameter values. 

4 Results 

4.1 Cortical response to COLA 

We first sought to quantify precisely the spatiotemporal dynamics of cortical non-muscle myosin II (NMY-2) that 
arises in response to ablating the cortex along a line. For this we used spinning disc microscopy to image NMY-2 
fluorescence in combination with a UV laser ablation setup to sever the cortical layer along a line 001321. We 
recorded the spatiotemporal evolution of the surrounding NMY-2 in the cortex to following the resealing process. 

In case of the zebrafish embryo, at the stage of 60% - 70% epiboly, COLA was performed within the YSL 
actomyosin ring and along a line parallel to EVL. In case of the C. elegans one cell embryo, the cortex was ablated 
just before the onset of cortical flow in the anterior half of the embryo and in a direction along AP axis. In both 
systems the area surrounding the cut was imaged until the hole was no longer visible due to turnover and regrowth 
(Fig.[2A and B right upper). 

COLA severs all connections within the cortex along the cut line and sets tension in the direction orthogonal to the 
cut line to zero. This results in tension gradients that drive an outward movement of the adjacent cortex away from 
the cut line (Fig.[T]A and B right). Notably, despite significant differences in cortical structure and dynamics, both 
zebrafish and C. elegans share large similarities in the overall response to COLA and the spatiotemporal dynamics 
of cut opening and resealing. In both systems the outward movement of the cut boundary lasts for several seconds 
and turns the cut line into an approximately elliptically shaped opening. Furthermore, the adjacent cortex moves 
outward with a velocity that decays over time, followed by cortical regrowth within the cleared region until no 
visible mark of the COLA procedure remains. 

To extract the characteristics of the response of the cortex to COLA and the characteristics of the resealing 
process, we analyzed the spatiotemporal dynamics of cortical NMY-2 after COLA in several ways. We first de¬ 
termined the outward velocity of cortical NMY-2 adjacent to the cut line by Particle Image Velocimetry (PIV) at 
a time just following the cut, as indicated by the arrows in Fig. [T] A and B right lower panels. Remarkably, the 
velocity profile along the direction perpendicular to the cut line was not uniform, but rather decayed over a char¬ 
acteristic distance away from the cutline as shown in Fig.^C and F. This spatial decay entails information about 
the characteristic distance A over which mechanical stresses communicate in the 2D cortical layer. Second, we 
quantified the temporal evolution of the cut opening. For this, we determined the extent of the opening generated 
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by COLA, by fitting an ellipse to this opening and determining the time evolution the minor radius of this fitted 
ellipse as a measure of the width of the cut opening (Fig.j^A). Notably, the width of the cut opening increased with 
time and reached maximum after approximately 3-5 seconds (Fig. |^D and G). The minor radius of the opening 
grows with a characteristic time governed by processes of stress relaxation in the actomyosin cortical layer. Third, 
we analyzed NMY-2 regrowth within the opening, by quantifying the average fluorescence intensity in a box that 
is placed at the center of the COLA opening (white broken line in Fig. A). Following an initial drop in intensity 
due to ablation and outward movement, NMY-2 intensity gradually recovered over approximately 30 seconds (Fig. 
I^E and H). Similar time scales have been observed in FRAP measurements that characterize NMY-2 turnover ||9|- 
These observations lend credence to the assumption that the cortex behaves as an active viscoelastic material. 
The shape evolution of the cut opening is determined largely by elastic properties of the cortex as well as active 
tension provided by NMY-2. The velocity decay away from the cut line is largely determined by the decay length 
of tension in the layer, and thus by 7 and A . Finally, re-association of myosin at the cut site is determined by 
actomyosin turnover. In what follows, we will compare the dynamics of opening and regrowth, as determined in 
our experiments by characterizing the spatial decay of the velocity held and the time evolution of the cut opening 
width and the myosin levels at the center of the hole, with theory. 




20 pm 




Figure 1; Probing cortical tension in the actomyosin ring of gastrulating zebrafish (A) and in the C. elegans 
(B)actomyosin cortex of single cell embryos by COLA. A Left, schematic of COLA (scissors) in the actomyosin 
ring of zebrahsh. COLA is performed along a 20 pm line (red line) at the stage of 65% epiboly. Red arrow-heads, 
direction of cortical recoil after ablation. Right, upper images show a time series of cortical NMY-2-GFP following 
ablation, lower images show the corresponding velocity helds as determined by PIV. B Left, schematic of COLA 
on the actomyosin cortex the C. elegans zygote. COLA is performed along a 10 pm line in parallel to AP axis of 
the embryo. Red arrow-heads, direction of cortical recoil after ablation. Right, upper images show a time series 
of cortical NMY-2-GFP following ablation, lower images show the corresponding velocity helds as determined 
by PIV. 
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Figure 2: COLA response in zebrafish and C. elegans : A and B, Illustration of the analysis of the cortical 
response to COLA. A, the time evolution of the cut boundary opening was characterized by fitting an ellipse to 
the cut boundary and measuring the minor axis of the fitted ellipse, white solid line. To analyze the regrowth of 
the cortical NMY-2 after the cut we measured the average fluorescence intensity within the box indicated by the 
white broken line. The size of the box is 1.8 x 1.4 pm for zebrafish and 1.4 x 1.2 pm for C. elegans. B To obtain 
the velocity profile along the direction perpendicular to the cut line (red line), the x components of the velocity 
vectors inside the rectangular area surrounded by the orange broken line were averaged in the y direction (box size, 
zebrafish; 27 x 13 pm , C. elegans: 12 x 10 jim). C and F, velocity profile along the direction perpendicular to the 
cut line at the time just after the cut for zebrafish (C) and C. elegans (F). D and G, time evolution of the minor axis 
of the fitted ellipse to the cut boundary for zebrafish (C) and C. elegans (F). E and H, average concentration of 
cortical NMY-2 as a function of time after the cut opening for zebrafish (C) and C. elegans (F). Red, experimental 
results in zebrafish (C-E ) and C. elegans (E-H). Error bars represents standard error with n=15 for all cases. Blue 
lines, theory curves utilizing the obtained least square fit parameters, see main text. 


4.2 Physical description of the actomyosin cortex 

We next sought to calculate the cut response in a coarse grained physical description of the cortical layer. Consid¬ 
ering that the thickness of the cortex is small compared to the size of the cell we describe the actomyosin cortex as 
an active 2D viscoelastic compressible fluid. We introduce a scalar field h{x, y, t) which denotes the local height 
of the cortical layer in the z direction, with x and y denoting spatial cartesian coordinates within the plane of the 
cortical layer, and t denoting time. Considering a viscoelastic isotropic active fluid in the plane and integrating 
over the height of the cortex we write the following constitutive equation; 

(1 -f tmA) {(T^J - (jfj) (3) 

— T Vh^k'^k^ij j 
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where indices *, j, k = x,y. The dynamic variables are the 2D stress tensor cr, the 2D active stress tensor and 
the 2D velocity field v. The material properties are characterized by the 2D shear viscocity rj, 2D bulk viscocity 
77b and a characteristic Maxwell time of stress relaxation tm- Here, Dt = dt + Vidi denotes the material time 
derivative. Turnover of the film is captured by the dynamics of the height given by 


dh 


-di{vih) + 



(4) 


where the first term on the r.h.s. accounts for advection of the gel by cortical flow, and the second term describes 
turnover which relaxes to the steady state hight h^. Note that both the turnover time Ta and the relaxed height /iq 
depend on actin and myosin turnover. The force balance equation reads fTl l 8 l l9l l3^ [34l 


= 'yvj , (5) 

where inertial forces have been neglected and 7 is a friction coefficient that describes frictional interactions be¬ 
tween the cortex and it its surrounding cytosol and membrane. 

We consider an active gel that is incompressible in 3D, with constant density through the height of the cortical 
layer. In this case, both the active stress and the 2D viscosities 77 and 77 b of the cortex are proportional to cortex 
height h and are written as 


r]{h) 

Vh{h) 


h 
ho 
Voh 
ho ’ 
Vhoh 
ho 


(6) 

(7) 

( 8 ) 


Here, denotes the isotropic active stress generated through ATP consumption of myosin, positive for 
contraction and dependent on the change in chemical potential associated with ATP hydrolyosis A 77 . Furthermore, 
rjo and yoo denote the shear and bulk viscosities of the layer when /i = /iq- 3D incompressibility condition couples 
divergences in the 2D flow velocity field v to change the cortex height h according to Eq.|^ and we set rjo/rjoo = 3. 
Eqs.[3]|^complete the model. Eor nondimensionalization we choose the characteristic time of regrowth Ta and the 
COLA cut length Ac as the respective time and length scales. Model equations with dimensionless variables can 
be rewritten with dimensionless parameters, ai = ^, a 2 = = ^ 7 ^?^ (see Appendix). 

We next asked whether this description can reproduce the relaxation dynamics in response to COLA that 
we observed in our experiments. To this end we numerically solved non-dimensional versions of Eq. [3]|^(see 
Appendix) in a rectangular box of width L ( ~ 15Ac) in x and y, with periodic boundary conditions. As an initial 
condition we choose uniform height and stress fields that correspond to the unperturbed stationary solution with 
V = 0. To account for COLA, this homogenous initial condition is perturbed at t = 0 by setting height h to zero 
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within a thin rectangular strip of length Ac and width ^ O.lAc. We then computed the resulting spatiotemporal 
dynamics. Fig.j^A displays the time evolution of the resultant height and velocity helds. We observe that 1) the 
velocity held is not uniform through the cortex but decayed over a distance from the cut line; 2) the width of the 
boundary initially grows and reaches a maximum; and 3) h{x^ y, t) in the the cut region recovers to steady state 
values on long times. These observations are consistent with those that we have made in our experiments, which 
lends credence to our approach. 

Next we analyzed our results from theory in terms of the spatial prohle of the velocity in the direction per¬ 
pendicular to the cut line, the growth of the cut boundary, and the recovery of the cut region in a manner that is 
similar to how we analyzed the experimental data. Fig. |^B to E demonstrates that the essential features observed 
in the corresponding graphs from our experiments (Fig. |^C to H) are reproduced in our theory. We next asked 
how changes of physical parameters of the cortical layer impact on the spatial prohle of the velocity in the direc¬ 
tion perpendicular to the cut line, the growth of the cut boundary, and the recovery of cortex height. To this end 
we performed numerical simulations of COLA with different values of a's. We hnd that increasing A leads to a 
corresponding increase of the spatial decay length of the velocity prohle in a direction perpendicular to the cut 
line (Fig.|^B), but has little impact on the growth timecourse of cut boundary (Fig.j^E). Eurthermore, increasing 
tm results in a corresponding increase in the relaxation time to reach the maximal width of the opening (Eig.j^C), 
with little effect on the hnal size of the hole. Taken together, key aspects of the relaxation process after COLA are 
separately determined by the characteristic length and time, A and tm of the cortical layer. 


A 



Distance from the cut (/im) 


Time (s.) 


Time (s.) 


Time (s.) 


Figure 3: Numerical calculation of the spatiotemporal response of a 2D active viscoelastic fluid in response 
to COLA: (A) Representative calculated nondimensional height (grey) and velocity held (arrows) of the cortex at 
times 0, 0.02,0.04,0.06, 0.08, and 0.68 following COLA by setting the height to zero in the dark region around 
the red line in the leftmost image (see main text). (B), Effect of variation of ai on the spatial velocity prohle 
at 0.5s after cut. The x component of velocity vectors was averaged along the y direction and within the box 
indicated by the broken orange line in the second panel of (A). (C-E), Effect of variation of a 2 and as on the 
time evolution of the half-width of the cut opening, determined by fltting an ellipse to the opening boundary. Note 
that changing ai has a small effect on time evolution of the width of the opening (E). Unless otherwise specifled 
parameter values are ai = 2.0, a 2 = 0.25, as = 4.5. Eor (B-E) we have dimensionalized time by Ta = 25s and 
length by Ac = 25y,m. 
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4.3 Comparison of theory and experiment for determination of physical parameters 


Next we determine coarse grained physical parameters of the cell cortex both in C. elegans single cell embryo 
and zebrafish actomyosin ring, by comparing the experimentally determined cortical response to COLA in terms 
of regrowth of the cut region, spatial decay of the outward velocity profile, and the time evolution of the cut 
opening boundary to theory. We note that the height variable h{x, y, t) in the theory is proportional to the per- 
area-concentration of molecules that are uniformly distributed along the z direction. Since the thickness of the 
cortical layer is smaller than the focal depth of the confocal microscope (approximately Ipm), the fluorescence 
intensity of cortical NMY-2, 1{x, y, t) in our microscope images is proportional to the 2D concentration of NMY- 
2, integrated over the height of the cortical layer. Hence the local 2D concentration I{x, y, t) of NMY-2 in our 
microscope images is proportional to h{x, y, t). We determined the best fit non-dimensionalized model parameters 
ai = ^, q ;2 = ^) c «3 = within a nonlinear least-square fitting scheme OTl by iteratively calculating the 

theoretical responses that best fit the experimental profiles (see blue curves in Fig. [^C to H in comparison to the 
experimental profile given by red curves). In Fig. [^C and D, we report the physical parameters of the cell cortex 
(Ta, A, tm, from the best fit model parameters. We obtained Ta = 23.4 ± 1.63s, A = 67.6 ± 8.2fj,m, tm = 
5.6±1.2sand ^ = 25.9±6.4/xm Vs for zebrafish, andra = 24.2±1.36s, A = 14.4±1.5/rm, tm = 4.48±1.25s 
and = 25.42 ± AAgLiri^/s for C. elegans. Note that these values were obtained by fitting each individual 
experiment (n=15 for both systems), and we report the respective averages ± standard error. Note also that the 
hydrodynamic lengths of ~ 68fim for the zebrafish actomyosin ring and ^ 14 /xto for the C. elegans cortex as 
well as the turnover times of ~ 25s show a close agreement with previous investigations ElISl. This supports our 
method of extracting values of physical parameters by use of fitting the measured COLA response to that expected 
from theory. 

The good agreement between theory and experiment indicates that the shape evolution of the cut opening 
might entail a sufficient amount of information for accurately determining both A and tm. Since the spatial 
profile of the velocity field is governed by hydrodynamic length. A, the temporal change of the shape of the cut 
opening boundary is likely to be affected by A. On the other hand, the temporal decay of the outward velocity is 
characterized by the time scale of the stress relaxation, tm. Thus it is possible that the shape evolution of the cut 
boundary is to a large extent governed by these two physical parameters. In agreement with this statement, the 
shapes of the cut opening boundary are distinct between the zebrafish actomyosin ring and C. elegans embryo, 
which might reflect the respective differences in physical parameters. 

Therefore we asked if it is possible to determine the hydrodynamic length A and the time scale of stress 
relaxation tm from the shape evolution of the cut opening boundary. To this end we determined the a’s by 
comparing the spatiotemporal development of the cut opening boundary shapes in experiment and theory. We 
detected the shape of the cut opening boundary observed in experiment to evaluate the difference with the shape 
from the theoretical prediction at the given parameter values. We computed the minimum distance from ith 
edge point detected in the experiment to the edge points obtained from the theory at the nth time frame, d™{n). 
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(see Materials and Methods for details). We then evaluated a merit function, D = where 

N is the number frames analysed (a useful choice for N was between 8 and 12 for a zebrahsh, and between 

4 and 6 for C. elegans). D was minimized iteratively to hnd the best ht values of a’s and Ta. Fig|^A and B 

shows the shape of the cut opening boundary from the theory with the best ht parameter values with respect to 
the minimization of D, successfully reproducing the shape evolution of the cut opening boundary observed in 
experiment. Consequently we obtained non-dimensionalized model parameters ai, a 2 , 0 : 3 , as well as the physical 
parameters Ta = 23.4 ± 1.66s, A = 80.2 ± 6 . 8 /j,to,tm = 4.95 ± 1.0s and = 26.89 ± 5.2/im^/s for 

zebrahsh and = 24.2 ± 1.33s, A = 13.6 ± 1.3pm, tm = 5.37 ± 1.2s and = 24.65 ± for 

C. elegans. Again, these values were obtained by htting each individual experiment (n=15 for zebrahsh and 
n= 10 for C. elegans, in the latter not all of our experimental datasets converged in the htting procedure), and 
we report the respective averages ± standard error. The values of the parameters of the actomyosin network 
are in close agreement between this method of determination, and the method used before, see Fig. [^C and D. 
We conclude that the shape evolution of the cut opening boundary entails a sufficient amount of information to 
determine the entire set of non-dimensional physical parameters, and provides a second means of determining 
physical parameters of the actomyosin cortical layer. 

5 Discussion 

While COLA is a powerful tool to characterize the tension in the actomyosin cortex, obtaining the physical pa¬ 
rameters to describe the cortical mechanics relies on an appropriate analytical method. Here we present a strategy 
to determine the physical parameters of actomyosin cortex from the cortical response to COLA. The method relies 
on the htting of the COLA response to the theoretical prediction from the 2D active viscoelastic huid model. We 
show that the best ht value of the hydrodynamic length, A’s are in close agreement with previous estimates both 
for zebrahsh actomyosin ring and C. elegans single cell embryo 121111 . Notably we can determine a whole set 
of parameters in a single experiment, in contrast the previous estimates that require the ensemble averaging for 
the how prohle, and myosin distribution. As shown in Fig. the standard error of the best ht parameters are 
small in both methods, signifying the overall robustness of the approach. In addition, our method does not require 
any assumptions for boundary conditions in the how and myosin prohles. Taken together, our method allows the 
precise determination of the physical parameters. 

Importantly, our method also allows to determine the characteristic time of stress relaxation, tm, which governs 
the time scale between elastic and viscous regime, tm sets the timescale for the large-scale movement of the 
actomyosin cortex, and morphogenetic deformations of the cortex that are driven by actives stresses in the layer 
generally occurs on time scales larger than tm. Recently active microrheology has been performed to measure 
the storage and loss moduli by manipulating a bead injected inside the cell 051 . This method allows the precise 
determination of the characteristic time of stress relaxation of the cytoplasm. However a typical diameter of bead 
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Figure 4: Determination of physical parameters by comparing the COLA responses observed in experiment 
to the computed response from the theory: A and B, Representative examples of fitting the cut opening bound¬ 
ary in response to COLA of the actomyosin ring in zebrafish during epiboly (A) and the cortex in single cell 
embryo of C. elegans (B). Red points, automatically detected boundaries of the cut opening. Distances between 
the computed and the detected boundaries were minimized to find the best fit parameter values. Cyan, the the¬ 
oretical boundaries that best fits the experimental ones. Images are one second apart. C and D, Comparison of 
physical parameter values between the two methods utilized. In red, the best fit parameter values were determined 
by comparing the experimentally determined regrowth of the cut region, spatial decay of the outward velocity 
profile, and the time evolution of the cut opening boundary to theoretical prediction (see Figure 2) . In blue, the 
best fit parameter values obtained by comparing the cut opening boundaries between experiment and theory, see 
(A) and (B). Note that both fitting procedures converge to similar values. Error bars are the standard error of the 
mean with n = 15 for zebrafish and 10 for C. elegans. 


used in the active microrheology measurement, ^ 500nm is bigger than the thickness of the cortical layer, one 
cannot distinguish the cortical and cytoplasmic viscoelastic properties in the measurement. In contrast, our method 
allows us to directly determine tm of the cortex. 

In summary, we have developed a robust and accurate method with broad applicability to determine the phys¬ 
ical parameters of cortical mechanics from COLA experiment in conjunction with the coarse grained mechanical 
theory. It provides us with the large-scale and biologically relevant parameters in terms of morphogenetic me¬ 
chanics. Given the simple description of the cortex used in the paper, we suggest that our method can be applied 
for complex multicellular systems such as epithelial tissues to determine the physical parameters that describes 
the tissue mechanics. 
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6 Appendix: Non-dimensional equations 


We rescaled the time and the spatial coordinate by setting the cut length (Ac) as length scale and as time scale 
of the system. Eqs 3]|5 can be rewritten as. 


^ = -a^dj (hd^atj) + 1- h 


(9) 

( 10 ) 

( 11 ) 


where A = represents the hydrodynamic length, ai = ^, a 2 = ^ and 03 = are three independent, 

dimensionless parameters of the model. 
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